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ABSTRACT 

Aims. Recent work has shown that major mergers of disc galaxies can only account for ~20% of the growth of the galaxy red sequence 
between z = 1 and z = 0. Our goal here is to provide merger frequencies that encompass both major and minor mergers, derived from 
close pair statistics. We aim to show that reliable close pair statistics can be derived from galaxy catalogues with mixed spectroscopic 
and photometric redshifts. 

Methods. We use B-band luminosity- and mass-limited samples from a Spteer/IRAC-selected catalogue of GOODS-S. We present 
a new methodology for computing the number of close companions, N c , when spectroscopic redshift information is partial. The 
methodology extends the one used in spectroscopic surveys to make use of photometric redshift information. We select as close 
companions those galaxies separated by 6/T 1 kpc < r v < 21h kpc in the sky plane and with a difference Av < 500 km s _I in redshift 
Q . space. 

Results. We provide A c for four different B-band-selected samples. It increases with luminosity, in good agreement with previous 
^ ' estimations from spectroscopic surveys. The evolution of N c with redshift is faster in more luminous samples. We provide A c of 

M* > 10 10 M Q galaxies, finding that the number including minor companions (A'" 1 , mass ratio yu > 1/10) is roughly two times the 
number of major companions alone (A,! 4 , mass ratio \i > 1/3) in the range 0.2 < z < 1.1. We compare the major merger rate derived 
by close pairs with the one computed by morphological criteria, finding that both approaches provide similar merger rates for field 
galaxies when the progenitor bias is taken into account. Finally, we estimate that the total (major+minor) merger rate is ~1.7 times 
the major merger rate. 

Conclusions. Only 30% to 50% of the M+ > 10 10 M Q early-type (E/SO/Sa) galaxies that appear between z = 1 and z = may have 
undergone a major or a minor merger. Half of the red sequence growth since z = 1 is therefore unrelated to mergers. 

£SJ ■ Key words, galaxies: evolution, galaxies: interactions, galaxies: structure 
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O ! 1- Introduction l2003b IGonzalez-Garcfa & Balcellsl l2005k iNaab et all l2006bl) . 

Observationally, merger remnan ts in the local Un i verse c an 
TT , , , , , . evolve into elliptical galaxies (Rothbere & Joseph 2006a b). 

How important are galaxy mergers on the mass assembly history In additiorii th e size |) addi et a lJ l2005h i Truiillo et al. 2007; 



(N 
> 



^ ; of red sequence (i.e., passive early-type) galaxies? This is one of gt ^ y£m dgr W el_etalJ S anTVelocity- 

~. ^ . the most challenging questions in galaxy evolution studies, one ~T • Tr$~^~^e~7~ t TTTTt7w* w e 

?H,. . ,t T i i • t,^ dispersion (Cenarro & Truullo 2009) evolution or massive 

that is motivated by both theory and observations. In the former, ^ redshiftj and ^ i u minosity density eV ol ut i on f 

popular hierarchical A-CDM models predict that the more ■ ... , jt, „ \ , |h n ^/iL WHC m 

F F . ,, ,. red sequence galaxies since z ~ 1 (Bell et al. 2004; Faber et al. 

massive dark matter haloes, inhabited by red-sequence galaxies, h m 4 , . ■ , .. , ■ , . 

_ , . . • ,, 2007) rule out the passive-evolution hypothesis and suggest 

are the final stage of successive mergers of less massive haloes. , . . • , c Z- 

, , e , . . , , 6 . . galaxy mergers are an important process in galaxy formation 

However, the behaviour of the baryonic component is still un- , „ 

_ , F , , , , and growth, 
clear. Only with many ad-hoc ingredients can the latest models, 

which include radiative cooling, star formation, and AGN and The merger fraction, / m , defined as the ratio between the 

supe rnova feedback, reproduce the observational trends better number of merger events in a sample and the total number of 

(see | Bower et alJfeOOd iDe Lucia & Blaizot 2007^ IStewart et alJ sources in the same sample, is a useful observational quantity 

2009bt iHopkins et al.l l2009cL and references therein). On the to explore the role of mergers in the growth of the red sequence 

other hand, A/-body simulations suggest that gas-rich merg- and thus to constraint cosmological models. Many studies have 

ers can produce intermediate-mass spheroida l systems (e.g., determined the merger fraction and its evolution with redshift 



Naab & Burkerj l2003t iBournaud et alJ 120051 : IHopkins et al.l up to z ~ 1, usually parametrized as f m (z) = / m (0)(l + zf , 
20081) . while dissipationless mergers can explain the more using different sample selections and methods, such as morpho- 



massive spheroids (e.g., IGonzalez-Garcfa & van Albadal logical criteria dConselice et al.ll2~003L I2008L l2009t lLaverv et al.l 



2 



C. Lopez-Sanjuan, et al.: Spectro-photometric major and minor companions in GOODS-S 



2004 ICassata et alJl2005b ILotz et al Il2008at IBridge et al.ll200 



Kampczv k et alJ 120071: iJogee et alJ 2009; Lopez-Saniua n et al 
2009a; Heiderman et al. 200 9tlDarg et al.l2010l). kinematic close 
companions (IPatton et all l2000i 120021: iPatton & Atfieldl l2008b 
ILin et all 12004 l2008t iDe Propris et al.1 120051 120071), s patiall 



(IFazio et alJl2004). we were able of re- 
separation (see IPerez-Gonzalez et alJ 



Lm et al.l lzUU4t IZUUift IJJ e Fropris et al .l IZUIDL KUU/I). spatially 

close pairslLe Fevre et alJ l2000t iBundv et all 12004" 2009a; 
u>-;^rrJ a t oi I nnml \v „r+„n ana a t oi I lonml lUo^v. Q t oi Tonne- 



Bridge et al 1 12001 IKartaltepe etail 120071: IHsieh et all 12008: 
Bluck et al .1120091). or the correlation functioii lBell et alj r2006; 
MasjedietaLll006). In these studies the value of the merger in- 
dex m at redshift z < 1 va ries in the range m - 0-4. A-CDM 
models predict m ~ 2-3 ( Kolatt et all 1 1 999t Governato et al.l 



1999: lGotflober et al.ll200ltlFakhouri & Mall2008l) for dark mat 
ter haloes, while suggesting a weaker evo l ution ; m ~ 0-2, for 
the g alaxy merger fraction dBerrier et al.l 120061: IStewart et al.l 
l2009ah . 

In this paper we study the number of close companions (N c ~ 
2/m) in a 5/9/fzer/IRAC-selected catalogue of the GOODS-S area 
and explore its dependence on redshift, B-band luminosity, and 
stellar mass. A robust methodolog y for measuring N c in spec- 
troscopic surveys was developed by P atton et aD (l2000l) . and we 
adapt it in present paper to exploit all the available redshift in- 
formation in spectro-photometric catalogues, obtaining reliable 
values when compared with those from fully spectroscopic sam- 
ples. In addition, work with close companions gives us useful 
information about the galaxies involved in the mergers, such as 
their mass ratio. Thanks to that, and taking advantage of our 
new methodology, we report an estimation of the number of mi- 
nor companions (stellar mass ratio 1:10) and their relation with 
major companions ( mass ratio higher than 1:3) for inter medi- 
ate mass galaxies (see lLotz et al.ll2008at IJogee et al.ll2009l for a 
morphological determination of the minor merger fraction). 

The paper is organized as follows. In Sect. [2] we summa- 
rize the GOODS-S data set we used, and in Sect. [3]we develop 
the methodology for determining the number of companions in 
spectro-photometric catalogues. Then, in Sect. 14. II we study the 
number of close companions in B-band selected samples, while 
mass-selected samples are analysed in Sect. 14.21 We compare 
our inferred major merger rates with those from morphologi- 
cal criteria in Sect. [5] We discuss the implications of our results 
in Sect. [6] and in Sect. [7] we present our conclusions. We use 
H Q = 70 km s -1 Mpc~', Q M = 0.3, and Q A = 0.7 throughout. 
All magnitudes are Vega unless otherwise noted. 



2. GOODS-S catalogue 

We worked with the galaxy catalogue from the Great 
Observatories Origins Deep Survey South (GOOD S-Sf] field 
by the Spitzer Legacy Team (Giavali sco et al.l 12004). We used 
the Version 1.0 catalogue^ and reduced mosaics in the FA35W 
(B435), F6Q6W (V 6 06), F115W (i 115 ), and FZ50LP (z 850 ) 
HST/ACS bands. These catalogues were cross-correlated us- 
ing a 1.5" search radius with the GOODS-S IRAC selected 
sample in the Rainbow co smological databasjl published in 
Perez-Gonzalez et al. (2008, see also lPerez^G onzalez et al.l2005l 
and Barro et al., in prep.), which provided us with spectral 
energy distributions (SEDs) in the UV-to-MIR range, well- 
calibrated and with reliable photometric redshifts, stellar masses, 
star formation rates, and rest-frame absolute magnitudes. We 
worked with a Spitzer/IRAC selected catalogue in [3.4] and [4.5] 
filters to ensure completeness in stellar mass. Although the PSF 



of Spitzer/IRAC is 
solve sources with ~ 1' 
2008 for details about the deblending process of Spitzer/IRAC 
sources). We used this 1" separation to fix the minimum radius 
when searching for close companions (Sect. l3~TT >. 

We refer the reader to these papers for a more detailed de- 
scription of the data included in the SEDs and the analysis pro- 
cedure. Here, we briefly summarize the main characteristics of 
the data set. The Rainbow database contains consistent aperture 
photometry in severa l UV, optical, NIR, and MIR bands with the 
method described in IPerez-Gonzalez et al (2008). The UV-to- 
MIR SEDs were built for 4927 IRAC sources in the GOODS-S 
region down to a 75% completeness magnitude [3.6] =23.5 mag 
(AB). These SEDs were fitted to stellar population and dust 
emission models to obtain estimates of the photometric redshift 
(Zphot), the stellar mass (M+), and the rest-frame B-band absolute 
magnitude (M#). 

Rest-frame absolute B-band magnitudes were estimated for 
each source by convolving the templates fitting the SED with 
the transmission curve of a typical Bessel-B filter, taking the 
redshift of each source into account. This procedure provided 
accurate interpolated B-band magnitudes including a robustly 
estimated ^-correction. Stellar masses were estimated using the 
expon ential star formation PEGASE01 models with a Salpeter 
(1955) IM F, and var ious ages, metallicities, and dust contents 
(ICalzetti et al.l 2000) were included. The typical uncertainties 
in the stellar masses are a factor of ~2 , which is typical of 
most stellar popula tion studies (see, e.g.. iPapovich et all [2006: 
iFontana et al]|2006l) . 

In the catalogue, ~ 40% of the sources have spectroscopic 
redshift (z spe c), and we rely on z p h t for the other ~ 60%. Because 
of this, we refer to our catalogue as spectro-photometric here- 
after. The median accuracy of the photometric redshifts at z < 
1.5 is |z spec - z D hotl/(l + Zsoec) = 0-04, with a fra ction <5% of 
catastrophic outliers (IPerez-Gonzalez et al. 1 120081 Fig. B2). In 
the present paper we use cr ZphM = 0-5A I + z pho t) as z phot error, 
where os„ is the standard deviation in the distribution of the 
variable 6 Z = (z p hot - z sp ec)/(l +z P hot), which is described by a 
Gaussian well with mean ux ~ and standard deviation as ; (see 
lLopez-Saniuan et al.l2009al for details). We take <x^ = 0.043 for 
z < 0.9 sources and os_ = 0.05 forz > 0.9 sources. 

Finally, we remove those sources in the catalogue within 
Az = 0.01 of the centre of the most prominent large-scale struc- 
ture (LSS) in GOODS-S located at z = 0.735 (iRavikumar et al.1 



1 http://www.stsci.edu/science/goods/ 

2 http://archive.stsci.edu/prepds/goods/ 
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http://guaix.fis.ucm.es/~pgperez/Proyectos/ucmcsdatabase.en.html 



120071: iRawat et al.ll2008t Lopez-Saniuan et al.ll2010h . This is be 
cause the relative velocity of two galaxies located in a cluster 
is representative not of the dynamical state of the pair, but of 
the cluster potential, and we cannot apply the Sect. [3]methodol- 
ogy. For that reason, the present results are mainly refer to field 
galaxies. 



3. Methodology 

In this sec t ion w e recall the methodology developed by 
Patto n et al. (2000), which has b een u s ed extensively on spec- 
troscopic samples (IPatton et al.1 120001, l200l IPatton ^ Atfieldl 
2008 : ILin et al.112004 l2008HDe Propris et al.ll2005l 120071 120 1 0b 
Ide Ravel et al.ll2009l) . Then, we extend that methodology to use 
all the available information in spectro-photometric samples, 
paying attention to different bias as luminosity/mass (Sect. 13. U 
and spectroscopic (Sect. 14. 1 . It completeness of the samples, the 
border effects in redshift space and images limits (Sect. [33}, and 
the treatment of multiple systems (Sect. l3~4l i. 
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3.1. Close pair statistics in spectroscopic samples 

The linear distance between two sources can be obtained from 
their projected separation, r p = 6dA(zd, and their rest-frame rela- 
tive velocity along the line of sight, Av = c \zi - Z;|/(l +Zi), where 
Zi and zi are the redshift of the primary (more luminous/massive 
galaxy in the pair) and secondary galaxy, respectively; 9 is the 
angular separation, in arcsec, of the two galaxies on the sky 
plane; and aU(z) is the angular scale, in kpc/arcsec, at redshift 



mm 
P 



< r r 



< r™ x and 



Z. Two galaxies are defined as a close pair if 
Av < Av max . The lower limit in r p is imposed to avoid identi- 
fying bright star-forming regions of the primary galaxy as close 
companions. Common limits are r™ n = 5h~ l kpc, r max = 20/i _l 

kpc, and Av max = 500 km s -1 . With these cons traints 50%- 
70% of the selected close pa i rs will finally m erge (Patton Tt alJ 
2000; iPatton & Atfieldl 120081: iLin et all 1200 4: Bell et al. 2006). 
We used the same limit in velocity but slightly different r p lim- 
its: r™ 11 = 6/T 1 kpc and r™ ax = 21/i _1 kpc. First, the minimum 
distance for which we are able to resolve two separate sources 
in our catalogue is 1" (Sect. [2]), which corresponds to ~ 8.5 kpc 
(6h~ l kpc) in our cosmology at the minimum of the function 
(Ja(z). This confusion limit fixes the value of r™ n . Second, we 
imposed r™ ax - r™" = 15/j" 1 kpc. This condition comes from 
the study of B ell et all ([2006). They find that the merger fraction 
is proportional to the radial range under study, so we kept the 
15/T 1 kpc range used widely in the literature. 

To compute close pairs we defined a primary and a secondary 
sample. The primary sample contains the more luminous source 
of the pair, and we looked for those galaxies in the secondary 
sample that fulfil the close pair criterion for each galaxy of the 
primary sample. 

If we work with B-band luminosity-selected samples, the 
primary sample comprises the galaxies in the catalogue with 



M, 



B,up 



< M B < 



M sel 



while the secondary comprises Mg up 



Mb < M S g \ galaxies, where Ms, U p is an upper limit in luminos- 
ity to avoid the different clusterin g properties of the most lumi- 
nous galaxies (Pat ton et al.ll2000l) . In every case, M s ^ 2 > M s £\. 
If we work with mass-selected samples, primary sample com- 
prises M s *\ > M+ sources, and secondary comprises M^ el 2 > M* 

f sel 



< M^\ . With the previous definitions 
the number of companions (^V c ) per primary galaxy is 



sources. In this case, M* 2 



1 N ' 

-7k 



(1) 



where N\ is the number of sources in the primary sample, and 
Nl the number of galaxies of the secondary sample that fulfil the 
close pair criterion for the primary galaxy i. Equation (Q]i is valid 
for volume-limited samples, but we work with lumino sity/mass- 
limited samples. To avoid incompleteness effects, Patto n et alJ 
(2000) define the function 



Sn(z) = 



(2) 



where Msximiz) is the limiting magnitude of the catalogue at red- 
shift z, Mg l = Mj| [Mg el 2 ] is the selection magnitude of the pri- 
mary [secondary] sample, and 0(Mb,z) is the luminosity func- 
tion in the Z?-band at redshift z. The definitio n of the function 
S jv(z) for a mass-limited sample is similar (see lRvan et all 2008. 
for details), and we take Sn(z) = 1 when Mi el < Msximiz), or 



Mt el > Mj,\i m (z). The hmitin g Mb magnitude was determined in 
Lopez-Sanjuan et al.j (l2009bl) : 



MbmM) = -13.78 - 12.66z + 11.18z 2 - 3.74z 3 , 



(3) 



defined as the third quartile in Mb distribution at each redshift 
dPerez-Gonzalez et a l. 2008). With this definition, the catalogue 
is complete for galaxies brighter th an Mb _ - 19.5 up to z ~ 1.3 . 
The limiting mass was defined in iPerez-Gonzalez et al .1 (I2008h 
as the 75% completeness of the catalogue for passively evolving 
galaxies, and can be parametrized in the range of interest as 



M*, lim (z) = 9.47z 1/8 M< 



(4) 



With this definition the catalogue is complete for galaxies more 
massive than M* = 6 x 10 9 M up to z ~ 1.3. 

We parametrize the luminosity/mass function in Eq. (f2]i with 
a Schechter function: 



<D(M,z) = Xiln(10)<f(z) 
x[10 X2 ' M * (z)_M) ] 1 +a< - z) exp[- 1 X2(M * (z)_M) ] 



(5) 



where M = Mb or M+, log 1() (0*(z)/Mpc 3 ) = O + j(z - 0.5), 
M*(z) = M* + S(z - 0.5), a(z) = a + if/(z - 0.5), and xi and 
X2 are constants to obtain the right normalization of the lumi- 
nosity/mass fun ction. We ob t ain th e B-band luminosity function 
parameters from Faber et al. (2007), while mass function param- 
eters from lPerez-Gonzalez et al.l d2008l) . For clarity we summa- 
rize all these parameters in Table Q] 

Finally, and following Patton et al. (2000), the number of 
companions normalized to a volume-limited sample is 



N r = 



^[■MzdZj-Mzj)" 1 ] 



Zf 1 S N (zd 



(6) 



where the index I covers all the close companions of the primary 
galaxy i. 

3.2. Close pair statistics in spectro - photometric samples 

The main problem in close pair studies with photometric sam- 
ples is to constrain the redshift space condition. For example, 
the Av < 500 km s~' condition at z,- ~ 0.7 implies \zi - z,| = 
Az < 0.0045. This condition is ~ 15 times less than the typical 
z p hot error at that redshift in our catalogue, cr Zj ~ 0.07. When 
one or both galaxies in a close spatial pair have photometric red- 
shift, we therefore cannot apply the methodology in Sect. 13.11 
To date a few works have used photometric catal ogues to de- 
termine pair statistics: Kartaltepe et alJ d2007l) and iBundv et al.l 
(l2009al) tackle the problem using a projection correction calcu- 
lated in random samples on the plan e of the sky, but ke eping the 
redshift information of the sources. iHsieh et all ((2008) redefine 
the velocity criterion to Az = 2.5cr 7 . and apply a conventional 
projection correction. Finally, iRvan et ail d2008l) also redefine 
the redshift criterion to Az = 2cr, , but no projection correction 
is applied. In this paper we present a new approximation to de- 
termine close pa irs in spectro - phot ometric samples based on the 
methodology of IPatton et al.l d2000l) and other previous photo- 
metric works. 

We use the following procedure to define a close pair system: 
first we search for close spatial companions of a primary galaxy, 
with redshift z\ and uncertainty cr zi , assuming that the galaxy is 
located at z\ - 2cr Zl . This defines the maximum 6 possible for 
a given rp 1 ™ in the first instance. If we find a secondary galaxy 
with redshift zi and uncertainty <x-, in the range r p < r™ ax and 
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Table 1. B-band luminosity and mass function parameters 



Sample selection 




x 2 


K 


6 


00 


r 


«o 





M B 


0.4 


0.4 


-21.07 


-1.23 


-2.46 


-0.12 


-1.30 





Mi, 


1 


-1 


11.23 


0.13 


-2.72 


-0.56 


-1.22 


-0.041 



NOTE. See Eq. (T5]l for details about the meaning of each parameter. 



with a given luminosity/mass with respect to the primary galaxy, 
then we study both galaxies in redshift space. For convenience, 
we assume below that every primary galaxy has, at most, one 
close companion in the secondary sample. In this case, our two 
galaxies could be a close pair in the redshift range 



[z~,Z + ] = [Zl - 2<T Zl ,Zl + 2cr z J n [Z 2 ~ 2<T Z2 ,Z 2 + 2cT z ,]. 



Because of variation in the range [z~,z + ] of the function cIa{z), 
a sky pair at z\ - 2cr Zl might not be a pair at z\ + 2cr Zl . We thus 
impose the condition r™ ln < r p < r™ ax at all z e [z~,z + ], and 
redefine this redshift interval if the sky pair condition is not sat- 
isfied at every redshift. After this, our two galaxies define the 
close pair system / in the redshift interval [zj,zj], where the in- 
dex j covers all the close pair systems in the sample. 

The next step is to define the number of companions asso- 
ciated at each close pair system j. For this, we suppose in the 
following that a galaxy i in whatever sample is described in red- 
shift space by a probability distribution P, (z,- 1 77,), where Zi is the 
source's redshift and 77, are the parameters that define the dis- 
tribution. If the source i has a photometric redshift, we assume 
that 



Pi (Zi I 7<) = p G (Zi I Zphoy. o", phoL .) = 



1 



V27T, 



■ exp< 



Trtr, 



(Zi Zphot,/) 

2^2 



while if the source has a spectroscopic redshift 

Pi (Zi I Hi) = P D (Zi I Zspec,,) = 5(Zi - Z spe c,i)> 



(8) 



(9) 



where 6(x) is delta's Dirac function. With this distribution we are 
able to statistically treat all the available information in z space 
and define the number of companions at redshift z\ in system j 
as 



Vj (Zl ) = Cj Pi(Zi I m ) jT " Pi(zi 1 772) dz 2 , 

where z\ e [zj,zj], the integration limits are 

z m =zi(l - Av max /c) - Av max / C , 
zi=zi(l+Av max /c) + Av max /c, 



(10) 



(11) 
(12) 



the subindex 1 [2] refers to the primary [secondary] galaxy in j 
system, and the constant Cj normalizes the function to the total 
number of galaxies in the interest range 



P\(z\ \rn)dzi + P2(z2\m)dz 2 



(13) 



Note that vi = if z\ < z , or z\ > zt- The function vi (Eq. IfTUl ) 

tells us how the number of close companions in the system j, N£, 
are distributed in redshift space. The integral in Eq. (TTOb spans 



those redshifts in which the secondary galaxy has Av < Av max 
for a given redshift of the primary galaxy. 

With previous definitions, the number of companions per pri- 
mary galaxy in the interval z r ,k = [z*,Z*+i] is 



(7) N t 



c.k 



Zjf Zk k+1 Vj(zi)dzi 

Zii Zk+i Pi(zi\m)<ki 



(14) 



where the index k spans the redshift intervals defined over the 
redshift range under study. If we integrate over the whole red- 
shift space, z r = [0, 00], Eq. (fT4l > becomes 



(15) 



where A^ c ' is analogous to N' c in Eq. ([TJ. 

The definition of function Vj in Eq. (TlOb is general, and we 
can find four different cases in our spectro-photometric samples: 

1 . Primary and secondary galaxies have z spe c- In this case Vj is 
V/(Zl) = 2X(J(Zl -Zspecl). (16) 

where z sp ec.i is the spectroscopic redshift of the more lumi- 
nous/massive galaxy in the close pair. With this definition, v ; 
is not zero only at z spe c,i- If all galaxies in the sample have 
z sp ec, Eq. (TOT l is equivalent to Eq. ([TJi. These systems provide 
N ] c -2 companions in Eq. (TOT l. 

2. The primary galaxy has z spe c and the secondary z p hot- 
Replacing the corresponding distributions of probability in 
Eq. (TlOl l we obtain 



v/(zi) = Cj6(z\ -z 



spec,l) ^ 



Pg(z2 I Zphot,2 , o- ZphoL2 ) dz 2 ■ (17) 



As in the previous case, Vj is not zero only at z spec ,i- These 
systems provide N% ~ 1 companions in Eq. ( T3"5l >. 
The primary galaxy has z p hot and the secondary z spe c- In this 
case Eq. (TT0T > becomes 



Vj(zi) = CjP c (zi |Zphot,l, cr, photl ). 

Function v ; (zi) is nonzero in the range 

Z S pec,2 - Av max /C Z sp ec ; 2+Av max /Cl 



1 + Av max /c 



1 - Av max /c 



(18) 



(19) 



This interval is imposed by the secondary galaxy, and spans 
those redshifts of the primary galaxy in which the secondary 
fulfils the condition Av < Av max . As in the previous case, 
these systems provide N[ ~ 1 companions in Eq. (fT5l >. 
Primary and secondary galaxies have Zphot- In this case the 
function v, is 

v j(Zl) = C/P G (ZlkphoU>0"z photi) 



r 



Pg(zi I Zphoa, o" Zphoa ) dz 2 . 



(20) 
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The main diffe rence between this approach and, for exam- 
ple, the one in iRvan et ail ((2008) is that we use the prob- 
ability distributions of the photometric redshifts to weight 
the number of companions in each system, N£, and to mini- 
mize projection effects. To illustrate how the weight process 
works, we show three examples in Fig.Q] In all these cases, 
the primary galaxy has z p hot,i = 0.6 and cr Zphotl = 0.03. In 
the panel (a), the secondary galaxy has z p hot,2 = 0.63 and 
cr Zphm2 = 0.03. With these values the function Vj is symmet- 
ric and not zero at [z~,Z^] = [0.57,0.66]. The integral of 
the function v, over the redshift space gives us the number 

of companions in the system, which is N£ = 1.64. In the 
panel (b) we increase the redshift of the secondary galaxy 
to Zphot,2 = 0.68 and keep its previous error. In this example 

the function v, is nonzero at [zT.zt] = [0.62,0.66], a nar- 

J j j 

rower range than in the previous case. Because of this, the 
number of companions in this system is only N£ = 0.46. 
Finally, in the panel (c) we keep the previous redshift of 
the secondary and increase its error to cr z = 0.06. In 
this case, the function Vj is not symmetric and not zero at 
[0.56,0.66], a similar range to the one in the 



fe7 



panel (a). However, the secondary galaxy is more extended 
in redshift space, and the number of companions is lower, 
Ni = 1.23. 

The last step involves normalizing our results to volume- 
limited samples, by applying the S # function (Eq. |2l). 



Vj(zi) = CjS N (zi)Pi(zi 



S N (Z2rP2(Z2\m)dZ2, (21) 



while the final number of companions is 



Zit^SxizdPiiZilTiddzi' 



(22) 



In ord er to estimat e the error of N c ,k we use the jackknife 
technique (Efr onll982l) . We compute partial standard deviations, 
6j, for each system j by taking the difference between the mea- 
sured N C £ and the same quantity with the jth pair removed for 



the sample, N* k , such that 5j = /V< 



c.k 



c,k 



For a sample with N 



systems, the variance is given by cr^ = [(A 7 - 1) J^jdA/N. 



3.3. Testing the method in a local, volume-limited sample 

We tested that our new methodology is able to statistically re- 
cover the number of companions in a spectroscopic survey from 
a photometric one. For th is, we study /V c in the Millennium 
Galaxy Catalogue (MGCQ, L iske et aljf2003l) . This survey com- 
prises 10095 galaxies with Bmgc < 2 over 37.5 deg 2 , 
with a spectrosc o pic co mpleteness of 96% dDriver etal . 2005). 
iDe Propris et alj d2007l) use the MGC to study the number of 
companions in a volume-limited sample (0.01 < z < 0.123, 
-21 < Mb — 5 log h < -18, A 7 = 3183 sources), obtain- 
ing /V C MGC = 0.035 ± 0.004 for rf* = and r™** = 20/T 1 
kpc (see also IDe Propris et al.l 120051) . We used this volume- 
limited sample in the present test. In addition, the MGC area 
had been observed by the Sloa n Digital Sky Survey (SDSS0, 
Adelman-McCarthv et al. 2006|), so every galaxy in the sample 



4 http://eso.org/~jliske/mgc/ 

5 http://sdss.org/ 
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Fig. 1. Probability distributions in redshift space of a primary 
galaxy (back solid line), a secondary galaxy (black dashed line), 
and the function vj/2 of the system (grey area). In all panels the 
primary galaxy has z p h t,i = 0.6 and (T Zpha] = 0.03. The sec- 
ondary galaxy has z p hot,2 = 0.63 and cr ZphoW = 0.03 in panel (a), 
Zphot,2 = 0.68 and cr. phoi2 = 0.03 in panel (b), and z P hot,2 = 0.68 
0.06 in panel (c). The number of close companions 



and cr 



2phol.2 



in each system, A^, is labelled in the panels. 



also has a photometric redshift. Comparing the z spe c's from MGC 
with the z p hot's from SDSS, we obtained cr, ta = 0.02. We take 
this uncertainty as representative of SDSS photometric redshifts. 

We defined / spec as the fraction of the sample's sources with 
spectroscopic redshift. The MGC sample has / spec = 1, while the 
SDSS sample has / spec = 0. To test our method at intermediate 
/ sp ec and for different cr. phoi , we assigned a z p hot to N(l - / spec ) 
random sources of the MGC sample, as drawn for a Gaussian 
distribution (Eq. (8j) with median z spe c and a given cr Zfm . Then 
we measured 



A/V c (/spec,Cr Zphol ) 



spec* O'jphot 



«) 



tfctt.O) 



(23) 
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Table 2. Number of companions in random samples in units of 
that in the spectroscopic MGC sample, AN C 



,/spec 


-phol 


= 0.01 


-phot 


= 0.02 


-■piioi 


= 0.03 





1.04 


±0.13 


1.24 


±0.12 


1.36 


±0.14 


0.2 


0.92 


±0.12 


1.06 


±0.13 


1.19 


±0.12 


0.4 


0.89 


±0.12 


0.93 


±0.12 


1.03 


±0.13 


0.6 


0.90 


±0.12 


0.93 


±0.12 


0.95 


±0.12 


0.8 


0.91 


±0.12 


0.96 


±0.12 


0.96 


±0.12 



where A^ c (/ spe c, o~z&a) is me number of companions in a random 
sample for a given / spec and cr Zphm , and N c (1,0) is the number 
of companions in the initial spectroscopic sample. We measured 
AAT C for /spec = 0, 0.2, 0.4, 0.6, 0.8, and o- Zpho , = 0.01,0.02,0.03. 
For each parameter combination, we repeated the process in ten 
different random catalogues and averaged the values. We sum- 
marize the results in Table [2] and show them in Fig. [2] We find 
that 

- When we only have spectroscopic information, we obtai n 
the same number of companions as De Pro pris et al.l (120071) . 
N c {\, 0) = 0.035 + 0.004. This implies that our methodology 
is equivalent to those used in spectroscopic samples when 
/spec = 1, as we hoped. 

- When we only have photometric information, our method 
recovers the initial number of companions (Fig. |2]i within 
error bars where the photometric redshift errors are small 
(cr Zphol = 0.01). On the other hand, the method overestimates 
the number of companions for cr Zphot = 0.02 and 0.03. This 
implies that we need small photometric redshift errors to 
avoid projection effects if only photometric information is 
available. 

- If we use the measured z p hot's of MGC sources from SDSS 
instead of the spectroscopic ones, the number of companions 
is higher than expected (A^ DSS = 0.041) and agrees with that 
from random samples with cr Zphoi = 0.02; that is, our random 
samples are representative of the observational ones, and we 
can use them to explore AA^ C at intermediate values of / spe c- 

- When our observational errors are not small enough, we 
need to increase / spec to obtain reliable N c values (Fig. 
That fixes the redshift of A^/ spec galaxies, and we use the 
photometric information to minimize projection effects. The 
higher the photometric redshift errors, the higher / spec must 
be to obtain reliable results. In this test we need / spec > 0.4 
for cr ZphM = 0.03 and / spec > 0.2 for cr ZplM = 0.02. 

Summarizing, our new methodology is able to provide re- 
liable N c values from spectro-photometric catalogues for either 
photometric redshift errors smaller than cr™ x or spectroscopic 

completeness higher than /™". For example, in our local test 
we find cr™ ax = 0.01, while /™" depends on the photometric 
redshift errors. These limits are not applicable to the GOODS- 
S catalogue because of the different galaxy number densities 
(e.g., the MGC vo lume-limited sample contains three clusters, 
iDriver et alJ[2005t while we remove the most prominent LSS in 
GOODS-S, see Sect. [2]) or the different search areas in the sky 
plane, which decreases by a factor of 15 from z = 0.1 to z = 0.7 
for the same physical area. Photometric redshifts errors in our 
GOODS- S cata logue are typically cr^ > 0.05, while / spec > 0.4 
(see Sect. l4.1.TI ). Our first step is therefore to compare our results 
at 0.2 < z < 1 .1 with those from spectroscopic samples to en- 




0.0 0.2 0.4 0.6 0.8 1.0 



/spec 

Fig. 2. AN C vs / spec for three different values of the photomet- 
ric redshift uncertainty, cr Zph01 = 0.03 (triangles; offset -0.05 for 
clarity), 0.02 (circles), and 0.01 (squares; offset +0.05 for clar- 
ity). We also show the value from MGC spectroscopic sample 
(/pec = 1; white pentagon and black solid line) and its 68% con- 
fidence interval (black dashed lines). The black star is the value 
obtained with the photometric redshifts of MGC galaxies from 
SDSS (/ spec = 0). 

sure that we can obtain reliable N c values from the GOODS-S 
catalogue. We perform this comparison in Sects. I3~4l and l4. II 

3.4. Number of companions as a function of r™ ax 

In this section we use the methodology developed in Sect. l3.2l to 
studying the variation in N c with r™ ax . This is a consistency test 

of our method because Nc oc rt~ y dPatton et al.l f2002: Lin et al.l 
2008; Ide Ravel et al.l2009l) . where y is the exponent in the corre- 
lation function, £(r) oc (rp/r) 7 . The y values in th e literature are 
consistent with y ~ 1.7 (e.g.. lLe Fevre et al .120051) . so we expect 
N c oc r p ' 3 . 

In this test we define two samples: one comprises all the 
galaxies in our catalogue with -22 < Mb < -19 and 0.2 < 
z < 1.1, named G p hot, while another comprises only the galax- 
ies of G p hot with z sp ec, named G spec - We take these two samples 
as primary and secondary, and we do not impose any difference 
in luminosity between both galaxies when searching for close 
companions. We vary r™ ax from 20h~ l kpc to 70h~ l kpc in 5/T 1 
kpc steps, and include the radius 21/i~' kpc. In all cases we use 
r™ n = 6/T 1 kpc. 

When we increase the radius of search, we start to find two or 
more secondary galaxies close to each primary galaxy. We treat 
these multiple systems with two different approaches: 

1. We study all the possible pairs as independent systems; that 
is, in a multiple system that comprises the galaxies A, B, and 
C; we study the subsystems A-B, A-C, and B-C. 

2. We only study the most massive pair, that is, the one with the 
lowest difference in luminosity/mass between the primary 
and the secondary galaxies. If this spatial pair is not a close 
pair in redshift space, we study the next more representative 
pair, and so on. 

We applied approach 1 to the G spec sample, while applying both 
to the G p hot sample, called G p hot,i and G p h t,2, respectively. We 
did this because approach 1 can fail in the G p hot sample owing 
to the high uncertainties in z p hot, making a projection correction 
necessary, but this is not the case if all the sources in the sample 
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Table 3. Number of companions as a function of r™ ax for ap- 
proaches 1 (independent pairs) and 2 (more representative pair) 
in the multiple-system treatment. 



20 30 40 50 60 
r max y h -\ kpc ] 



70 



Fig. 3. Number of companions vs r™ ax for galaxies with -22 < 
Mb < -19 and 0.2 < z < 1.1. White dots are obtained using ap- 
proach 1 in the multiple-system treatment (G p hot,i), while black 
dots are from approach 2 (G p hot,2, see text for details). White 
squares are from approach 1 over sources with spectroscopic 
redshift (G speCi i). The black solid line is the best power-law fit, 
N c oc r^, to the G S pec,i data, grey solid line to the G p hot,i data, 
and grey dashed line to the G p hot,2 data. The exponent of the fits 
is labelled in the figure. The black dashed line is the expected 
number of companions in the G p hot sample, derived from G spec ,i 
data. 

have z spec . Approach 2 is more conservative and we can miss, for 
example, a minor companion of a galaxy with a major one, but 
mitigate the contamination by projection effects. We summarize 
the results in both Tableland Fig. [3] Fitting a power-law func- 
tion to the data, N c oc r p , we obtained q = 1.65 in the G p hot,i 
case, q = 1.32 in the G p h t,2 case, and q — 1.21 in the G spec ,i 
case. We see that 

- The trends obtained in the G spec ,i and G p hot,2 cases agree with 
the expected q ~ 1.3; however, the value of q in the G p hot,i 
case is higher than expected, reflecting the increased contam- 
ination by projection effects when the search area increases. 

- For a given r p nax , one expects that N c oc n g , w here n g is 
the number of galaxies in the sample dLin et al . 2008); in 
our case, «g(G spec ) = 0.47n g (G p h o tX so we can predict the 
expected number of companions in the G p hot sample from 
the obtained in the G spe c,i (Fig.©. We see that the expected 
number is in excellent agreement with the G p hot,2 number of 
companions. This implies that our methodology is able to re- 
cover statistically, from the GOODS-S spectro-photometric 
sample, the same N c that in a fully spectroscopic one, as we 
hoped (see also Sects. [3~3l and l4. 1 .Tb . 

- The number of companions at rL™* < 25h kpc in the G p h t,i 
and G p hot,2 cases are the same. This implies that projection 
effects start to be important at r™" > 30/i~' kpc. 

In this section we have applied our new methodology to the 
GOODS-S catalogue with positive results, and in the following 
we use approach 2 to treat multiple systems. Note that 47% of 
the G p hot sources have z spec . We study how the spectroscopic 
completeness of the sample affect our results in Sect. I4.1.T1 (see 
also Sect[13l. 



3.5. Border effects in redshift and in the sky plane 

When we search for a primary source companion, we define a 
volume in the sky plane-redshift space. If the primary source is 



..max 

r p 

(/r 1 kpc) 




^.(0.2 < z < 1.1,- 


22 < M 


B <-19) 




G 


3hot,l 


Gphot,2 


^spec.l 


20 


0.049 


± 0.006 


0.049 ± 


0.006 


0.028 ± 


0.008 


21 


0.056 


± 0.006 


0.056 ± 


0.006 


0.030 ± 


0.008 


25 


0.081 


± 0.007 


0.077 ± 


0.007 


0.038 ± 


0.009 


30 


0.112 


± 0.008 


0.101 ± 


0.008 


0.048 ± 


0.010 


35 


0.148 


± 0.009 


0.128 ± 


0.008 


0.060 ± 


0.011 


40 


0.186 


± 0.009 


0.152 ± 


0.009 


0.065 ± 


0.012 


45 


0.217 


± 0.009 


0.169 + 


0.009 


0.075 ± 


0.013 


50 


0.248 


± 0.009 


0.190 + 


0.009 


0.088 ± 


0.013 


55 


0.290 


± 0.009 


0.219 + 


0.009 


0.098 ± 


0.014 


60 


0.326 


± 0.008 


0.231 ± 


0.009 


0.103 + 


0.014 


65 


0.367 


± 0.008 


0.251 ± 


0.009 


0.117 + 


0.015 


70 


0.399 


± 0.008 


0.257 ± 


0.009 


0.136 + 


0.016 



near the boundaries of the survey, a fraction of the search vol- 
ume lies outside of the effective volume of the survey. To ac- 
count for this, we use a correction factor ft,, the fraction of the 
area 7r[(r p nax ) 2 - (r™ 11 ) 2 ] around the primary galaxy that lies in 
the survey area. This factor is positive and depends on z because 
the projected distance in the sky is function of primary source 
redshift. We studied how this factor affects the previous section 
results, finding that border effects are representative (i.e., Icr dis- 
crepancy) only at r™* > 70/i~' kpc. 

We avoid the redshift incompleteness by including in the 
samples not only the sources inside the redshift range under 
study, but also those sources with zi +2<x,- > 0.2 and zi-lcri < 1 . 1 
(see Sects. [3~4l and l4.U . This implies that our samples comprise 
sources with z e [0.1, 1.3). 



4. Results 

4. 1 . Number of close companions in B-band 
luminosity-selected samples 

In this section we study the number of companions as a function 
of z and Mb- For this, we define three luminosity-selected sam- 
ples, -22 < M B < -20, -22 < M B < 19.5, and -22 < M B < 19, 
and study A^ c in two redshift ranges, named z r ,i = [0.2,0.65) 
and z r ,2 = [0.65, 1.1). We chose these redshift ranges to ensure 
good statistics in our study. In addition, we take the three previ- 
ous luminosity-selected samples as primary and secondary, and 
do not impose any limit in primary to secondary luminosity ra- 
tio. We apply the same luminosity selections and search param- 
eters to the MGC sample described in Sect. 13.31 to obtain con- 
sistent data points at z = 0.092. We summarize our results with 
those from previous spectroscopic studies in literature in Table|4] 
and show them in Fig. [4] We find that our results from spectro- 
photometric samples are in excellent agreement with those from 
fully spectroscopic ones when similar search parameters (i.e., 
luminosity selection and radius range) are applied. 

The dependence of A^ c on z can be parametrized with a 
power-law function, 



N c (z) = JV C (0)(1 + zf 



(24) 



The least-squares fits of Eq. (l2"4t to the data are summarized 
in Table [5] We find that the number of companions evolves 
faster for more luminous galaxies. The index m decreases from 
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Fig. 4. Number of companions vs redshift for different fi-band luminosity selections: panel (a) for -22 < Mb < -20 galaxies, panel 
(b) for -22 < M B < -19.5 galaxies, panel (c) for -22 < M B < -19 galaxies, and panel (d) for -22 < M e B (AB) < -19 galaxies 
(see text for details ). Symbols are the n umber of com panions from this work ( GOODS-S. whit e squares; MGC, white t riangles), 
Patto n et ail ([2000, dark grey squares), iPatton e t al. (2002, dark grey circles), iLin et al.l ([2004, light grey circles), and iLin et al] 
(2008, light grey triangles). Errors bars in z axis show the redshift range spans for each point. Black solid lines in all the panels are 
the best least-squares fit of a power-law function to the da t a. The parameters of the fit are in each panel. Black dashed line in panel 
(d) is the power-law parametrization provide bv lLin et"al I (120081) with m = 0.4. To facilitate comparison, the scales are the same in 
all the panels. 



Table 4. Number of companions in B-band luminosity-selected samples 



Reference 


z 


-22 < M 


5 <-20 


-22 < M B 


< -19.5 


-22 < M 


s < -19 


Present work (GOODS-S) 


0.425 


0.022 ± 


0.013 


0.029 ± 


0.010 


0.036 ± 


0.009 




0.875 


0.041 ± 


0.010 


0.057 ± 


0.009 


0.064 ± 


0.008 


Present work (MGC) 


0.092 


0.011 ± 


0.004 


0.024 ± 


0.005 


0.034 ± 


0.005 


Patton et al. (2000) (a) 


0.015 


0.008 ± 


0.002 


0.012 ± 


0.003 


0.019 ± 


0.004 


Patton et al. (2002) (b) 


0.297 


0.014 ± 


0.003 


0.023 ± 


0.006 


0.035 ± 


0.008 


Lin et al. (2004) (c) 


0.1 
0.6 
1.1 


0.049 ± 


0.013 


0.036 ± 


0.011 


0.024 ± 


0.010 


Lin et al. (2008) (c) 


0.077 
0.48 






0.036 ± 


0.005 


0.031 ± 


0.004 




0.85 


0.040 ± 0.005 











NOTES, (a) Original data are obtained with h = 1. (b) Original data are obtained with h = 1 and M B (z) = M B — z sampl e selec tion, (c) We use 
the Eq. (T24T> par ameters provide for -21 < M B (z) (AB) < -19 galaxies, where M e B (z) - M B - Qz, with Q = 1 in lLin et all (2004) and Q = 1.3 in 
ILin et al] (12008). The searc h radius is 10/r' < r p < 30/r 1 kpc, so we multiply their values by 3/4 to normalize their results at 15/r 1 kpc (from the 
correlation funtion study of lBell e t al. 2006). 



m = 2.5 + 0.1 for M B < -20 galaxies to m = 1.6 + 0.4 
for Mb < -19 galaxies. This trend is well described as m cc 
(-0.8 + 0.4) x Mf. On the other han d AT c (0) increases when lu- 
mino sity decreases, i n agreement with Patton et al. (2000, 2002); 
and lLin et al. (2004) results. These trends also agree with those 



predic ted by the cosmological model of Kho chfar & Burkertl 

(12001b . 

In the near universe (z ~ O.D. lPatton & Atfieldl (120081) study 
the dependence of Af c on r-band luminosity in an SDSS sam- 
ple. They find a nearly constant value over the range -22 < 
M r < -18 for r™ n = 5/1" 1 kpc and r™ ax = 20/1" 1 kpc close 
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companions, = 0.021 + 0.001, in contrast to our measured 
evolution with luminosity of N C (Q). However, they only account 
for major companions (AM, = 0.75), while we do not impose 
any luminosity constraint. This implies that we are sensitive to 
minor companions, more numerous than major (see Sect. 14. 21 , 
in the lower lumin osity sample, that le ad to an increase in N c . 
On the other hand, Ide Ravel et al.l (120091) measure the evolution 
of the major merger fr action (AM B = 1.5) in the VIMOS-VLT 
Deep Survey (VVDS0, Le Fevre et al .]|2005|) up to z ~ 1, finding 
that the index m is lower in more luminous samples, an opposite 
trend to what we find. This discrepancy can be explained again 
by the different definition of companion: we lose major compan- 
ions near the selection luminosity, but gain the minor ones of 
the more luminous galaxies. These two examples point out that 
similar definitions of companion are needed to compare results 
from different studies and surveys. 

iLin et all d2008) study the evolution of N c at < z < 1.2 
from ~ 35000 spectroscopic sources, finding A^ c (0) = 0.031 + 
0.004 and m - 0.4 ± 0.2, a low value that is in c ompa tible with 
those in the present work. However, ILin et al galaxies 
have -21 < M e B (z)(AB) < -19, where M e B (z) = M B - 1.3z, 
a different selection than ours. To explore whether this can 
be the origin of the d iscrepancy in m value, as sug gested by 
iKartaltepe et al.l d2007l) . we mimic the lLin et al.ld2008l) selection 
in the lite rature and in pre sent work as close as possible. In the 
latter, the lLin et al"1(l2008l) selection at z = 0.425 is M B < -19.5 
and is M B < -20 at z = 0.875. We summarize the data in Table|6] 
and show them in the panel (d) of Fig.|4] We obtain m = 0.7+0.4, 
a low value compatible with the ILin et al.l (120081) result: the se- 
lection of the sample is a key issue in determining m, and must be 
taken into acc ount when one com pares different works. Finally, 
the models of lBerrier et al.l (|2006) predict m — 0.4 - 1 .0, in good 
agreement with our result. 



4.1 .1 . Dependence on the spectroscopic completeness of 
the sample 

In previous section we show that our methodology provides reli- 
able N c values, compatible with those obtained in fully spectro- 
scopic samples, from GOODS-S spectro-photometric samples. 
As we show in Sect. l3.3l this can stem from either i) small photo- 
metric redshift errors or ii) enough spectroscopic completeness. 
The -22 < M B < -20 sample has / spec = 0.54, where / spec is 
the fraction of sources with z spec in the sample, the -22 < M B < 
-19.5 sample has / spec = 0.49, and the -22 < M B < -19 sam- 
ple has / sp ec = 0.42. This means that our methodology works at 
least for / spec > 0.4 samples. To check that photometric redshift 
errors in GOODS-S catalogue are small enough to skip spectro- 
scopic information, we repeat the study in the previous section 
but using z p hot for all the sources, although some have z spec (i.e., 
/spec = 0). Unfortunately, we obtain higher values ofN c than ex- 
pected. This implies (i) that we need / spec > 0.4 in the current 
GOODS-S samples to avoid an overestimation in A^. because of 
projection effects and (ii) that lower as, are needed in order to 
apply our methodology to / spec — catalogues up to z ~ 1 (e.g., 
the COSMOS0 survey, where cr Sz ~ O.Ol. lIlbert et all[2009h . 



6 http://www.oamp.fr/virmos/vvds.htm 

7 Cosmological Evolution Survey, IScoville et al.1 120071 
(http://cosmos.astro.caltech.edu/index.html I. 



Table 5. Fit parameters of the function N c (z) = N c (0)(l + z) m to 
the data 



Sample selection 


/Vc(0) 




m 


-22 <M B < -20 


0.008 ±0.001 


2.5 


±0.1 


-22 < M B < -19.5 


0.014 ± 0.002 


2.2 


±0.3 


-22 <M B < -19 


0.023 ± 0.003 


1.6 


±0.4 



Table 6. Number of companions for -21 < M e B (AB) < -19 
galaxies 



Reference 


z 


K 




Lin et al. (2004) 





0.022 ± 


0.010 


Lin et al. (2008) 





0.031 ± 


0.004 


Patton et al. (2000) 


0.015 


0.021 ± 


0.005 


Present work (MGC) 


0.092 


0.034 ± 


0.005 


Patton et al. (2002) 


0.297 


0.028 ± 


0.007 


Present work (GOODS-S) 


0.425 


0.029 ± 


0.013 


Present work (GOODS-S) 


0.875 


0.041 ± 


0.010 



4.2. Number of close companions in mass-selected samples 

In this section we study the number of major and minor com- 
panions of primary galaxies with log(M^\/M e ) = 10. As in 
previous sections, we define two redshift ranges, named z r ,\ = 
[0.2,0.65) and z r ,2 = [0.65, 1.1). In the following we denote the 
galaxy mass ratio as 



where M+ \ and M*^ are the stellar mass of the primary and the 
secondary galaxies in the pair, respectively. To explore the de- 
pendence of N c on p, we vary the mass ratio under study from 
p > 1 /2 to p > 1/10. We take major companions, denoted N™, 
as those with p > 1/3, while minor companions, denoted N"\ 
have p > 1/10. With these definitions, major companions are 
included in A 7 ™. For completeness, we take = p x 

hereafter, unless noted otherwise. All the mass-selected samples 
in this section have / spec > 0.4 (see Sect l4.1.j"l for details), vary- 
ing from / spe c = 0.65 for M* > 10 10 M galaxies to / spec = 0.41 
for M* > 10 9 M Q galaxies. 

We summarize our results in Tableland show them in Fig. [5] 
We find that (i) the number of close companions decreases with 
redshift for every fi. Combining all the fi values, this evolution 
is described by N c ( z) °c (1 +^) 2 - 1±Q - 5 well . This evolution agrees 
with the finding by Ide Ravel et al. I (120091) for the major merger 
(/J. > 1/4) fraction of M* > 10 10 M Q galaxies in VVDS-Deep 
survey, m = 2.04. And (ii), the number of close companions 
grows when fi decreases. Combining both redshift ranges, we 
obtain N c m = (1.7 ± 0.3) x N™; that is, the number of minor 
companions is roughly twice the number of major companions. 

Only a f ew wo rks have studied minor mergers statistically. 
Joge e et al.l (12009) report the minor merger (p > 1/10) frac- 
tion in Galaxy Evolution from Morphology and S EDs (GE MSEl 
iRix et al.l l2004) for M* > 4 x 10 10 M G galaxies dSalpeterll 19551 
IMF), and estimate tha t minor mergers a re three times ma- 
jor (p > 1/4) mergers. ILotz et al . (2008a) use G - M20 mor- 
phological indices to determine the minor merger (p > 1/9, 

8 http://www.mpia.de/GEMS/gems.htm 
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Table 7. Number of close companions of M^i > 10 10 M galax- 
ies as a function of mass ratio fi 



Fig. 5. Number of companions vs mass ratio fi for primary galax- 
ies with > 10 10 M Q . Open circles are forz = 0.875 galaxies, 
and grey squares for z = 0.425 galaxies. The solid lines are the 
best fit of a power-law function, N c oc //, to the data. The black 
dashed line is the best fit of a power-law function with a fixed 
exponent s = -0.6 to the z = 0.425 data (see text for details). 



lLotz et al.1 120091) fraction in All-Wav elength Extended Groth 
Strip International Survey (AEGIS0, iDavis et al.ll2007l). T heir 
values are in good agreement with the iJogee et all d2009l) re- 
sults, but they do not estimate the major merger fraction. In the 
local Universe. iDarg et aD (12010) estimate that minor merger s 
are twice major mergers in Galaxy Zoc[3 dLintott et al.l 12008): 
the latter is based on the visual classification of SPSS galaxies 
by internet users. Finally. IWoods & Gellerl (l2007h study the dif- 
ferent properties of major (Am, < 2) and minor {Am, > 2) close 
pairs in SDSS, which corresponds to fx ~ 1 /7. Unfortunately, 
they do not attempt to derive merger fractions, but the influence 
of clo se companions on galaxy properties (see also lEllison et al.l 
2008). Summarizing, literature values are consistent with N™ — 
2 - 3 x N™. The different methodologies (close pair vs mor- 
phology) and sample selections make quantitative comparisons 
difficult (see also Sect. 16. U . but the qualitative agreement is re- 
markable. 

The observed dependence of N c on \i is para metrized well 
as N c ( iu) oc [X s , as predicted by the simulations of iMaller et al.l 
(2006). Fitting a power-law function to the data we obtain s = 
-0.6 + 0.2 at z = 0.875, and s = -0.4 + 0.2 at z = 0.475. 
However, the observed evolution in s is not significant: if we im- 
pose s = -0.6 in the range z r .\, the observational values are also 
described well within the uncertainties. The highest difference 
between the model with fixed s = -0.6 and the observations 
occurs at /j > 1/10, i.e., for minor companions. This suggests 
that log(M*/M ) > 10 galaxies have accreted satellite galaxies 
between z ~ 0.9 and z ~ 0.5. Since most of these low-mass 
satellites are expected to be gas-rich, this accretion could ex- 
plain the resid ual star formation o bserved in early-types galax- 
ies at z ~ 0.6 (Kavirai et al. 2010). Studies with larger samples 
are needed to constraint the differential evolution with redshift, 
if any, of the number of major and minor companions. 



Secondary sample 


A< 


z = 


0.425 


z = 


0.875 


log(M*, 2 /M ) > 9.7 


1/2 


0.018 


±0.011 


0.029 


± 0.009 


log(M*, 2 /M Q ) > 9.5 


1/3 


0.028 


±0.014 


0.039 


±0.011 


log(M*, 2 /A/ ) > 9.4 


1/4 


0.031 


±0.015 


0.046 


±0.011 


logtM^/Me) > 9.3 


1/5 


0.032 


±0.015 


0.060 


±0.013 


log(M,, 2 /M e ) > 9.1 


1/8 


0.038 


±0.016 


0.075 


±0.014 


log(M*, 2 /M Q ) > 9 


1/10 


0.038 


±0.016 


0.079 


±0.014 



5. Major merger rate: close pairs vs morphological 
criteria 



The local study (z ~ 0.09) of lDePropris et al l (120071) shows 
that merger fractions by close pairs and by morphological crite- 
ria (i.e., taking highly distorted galaxies as major merger rem- 
nants) give similar merger rates when samples are compared 
carefully. However, in the range 0.2 < z < 1.2 both methods 
yield diffe rent merger rate s, some times by an order of magni- 
tude (e.g.. iLin et al.ll2004l) . In this section we compare the ma- 
jor m erger rate inferred by our clo se pair study, , with that 
from Lopez-Sanjuan et al. (2009b, L09 hereafter) by morpho- 
logical criteria, 5l^ ph . In L09 gas-rich major merger remnants 
are select ed as those galaxies with high values o f the asymmetry 
index (A. lAbraham et alJ[l996 : Conselicel l2003l) . L09 determine 
the morphological merger fraction of log(M*/M ) > 10 galax- 
ies in the same GOODS-S catalogue that we use in the present 
paper, and lead with three important sources of systematics: (i) 
they avoid morphological /f-corrections by measuring the asym- 
metries in the rest-frame B-band; (ii) the deal with the loss of in- 
formation with redshift (i.e., spatial resolution descent and cos- 
mological dimming) by artificially redshifting all the galaxies to 
a unique and representative redshift, Zd — 1 ; and (iii) they take 
the effect of observational errors in z and A i nto account by max- 
imum likelihood techniques developed in Lopez-Sanjua n et al.l 
j2008l). To obtain %." ph from the merger fraction, they assume 
Perez-Gonzalez et al I (120081) mass functions, as in the present 
paper , and a typical timescale of ~ 0.5 Gyr (Conselice 2006, 
[20091: lLotz et al.ll2008bl l2009h 

Following lLin et all d2004l) . we define the major merger rate 
by close pairs as 

M*) = ~C m P (z, M*) Nf{z, M±) r p -' r , 



(26) 



9 http://aegis.ucolick.org/ 

10 http://www.galaxyzoo.org 



where the factor 1/2 is for obtaining the number of merger 
systems from the number of close companions, the factor 4/3 
takes the lost c ompanions in the inner 6/T 1 kpc into account 
dBell et al.l [2006). C m accounts for the fraction of the observed 
systems that really merge in r pa i r , and p(z, M+) is the comoving 
number density of galaxies more massive than at redshift 
z. We take C m = 0.6 ± 0.1 dPatton et al.ll2000l: iBell et alJl2006l: 
ILin etal.ll2008l:lde Ravel et all 120091). 7 Da j r = .75 + 0.25 Gyr 
dKitzbichler & Whitd2008Hde Ravel et al.|2009t). and determine 
p(z, M +) with the mass functions from IPerez-Gonzalez et al.l 
(120081) . 

Although L09 determine ^R™ ph for M* > 10 10 M Q galax- 
ies, we cannot compare their merger rate with that from 
N™(z, 10'° Mp) b ecause of the progenitor bias (IBell et al.l 2006: 
lLotz et al]|2008al) . In morphological studies we are sensitive to 
the high distorted remnant phase of a major merger (jj > 1/3, 
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IConselicell2006t lLotz et aD 12009). while in pair studies we see 
the pre-merger stage; that is, the mass of the future remnant is the 
su m of the two galax ies' masses in the pair. This also was noted 
bv lGenel etal.1 (120091) . who compare the merger rate per progeni- 
tor (i.e., close pairs) and descendant (i.e., morphol ogical criteria) 
dark matter halo in the Millennium Simulation ( Springel et al. 
120051) . They find that these two merger rates are different quan- 
tities, so close pairs and morphological studies do not measure 
the same merger rate. We therefore compare the morphologi- 
cal merger rate of M* > 10 10 M G galaxies with that from 
obtained with ^(M^/Mq) = 9.7, log(Mf 2 /M G ) = 9.4, and 
jj > 1/3. With this definition, the merger of two galaxies with 
limiting masses is log(M*/M ) ~ 9.9, a mass that could in- 
crease o wing to the induced star formation in gas-rich merg- 
ers (e.g.jLin et al. UgOOl iLi et al.ll2008l iKnapen & James 2009; 
Robaina et al. 2009). Finally, we apply a r pa i r = 0.75 Gyr de- 
lay to the close pairs data points to mimic the redshifts at which 
the highly distorted remnants of these close pairs systems could 
be observed. The higher redshift point at z — 0.875 becomes 
Z = 0.725, while the point at z = 0.4 becomes z = 0.335. In sum- 
mary, both morphological and close pair studies now provide the 
rate of major merger remnants with M+ > 10 10 M Q (i.e., per de- 
scendant galaxy), and we can compare each of them. 

We summarize the final values of from this work and 

those of 9lJ^ ph from L09 in Table [8] and show them in Fig. [6] We 
see that both methodologies provide similar merger rates within 
the error bars in the range under study, as expected if we are ob- 
serving two different phases of the same physical process in the 
same field (i.e., we minimize the field-to-field variance between 
both measurements). This reconciles the morphological merger 
rates with those from pair statistics and lends credibility to the 
asymmetry index A as a major merger indicator. 

We parametrize the merger rate evolution with a power-law 
function, 



R m = Km(0)(1+z)". 

The fit to all the data in Table |8]yields 



(27) 



(28) 



^R M = (0.42 ± 0.07)(1 + zr s±u ' 3 x 10~ 4 Mpr 3 Gyr" 1 . 

These parameters are substantially different from those in L09, 
who find % M (0) = (0.29 ± 0.06) x 10~ 4 Mpr 3 Gyr 1 and n = 
3.5 + 0.4 from morphological data alone. 

In the previous discussion we do not consider that the close 
pair methodology is sensitive to gas-rich (wet + mixed) and 
spheroidal-spheroidal (dry) mergers, while the mor phological 
merger rate only refe rs to gas-rich major mergers dConselicel 
2006; Lot z et al.ll2009t Lopez-sanjuan et al., in prep.). However, 
at z > 0.3 th e fraction of dr y merge rs with respect to th e to- 
tal is ~ 10% dLin et al.l l2008. see also Ide Ravel et al.ll2009l) . in- 
creasing their importance at z < 0.2 and being ~ 25% of the 
total at z ~ 0.1 dLinet al.ll2008l This implies that % p ™ must 

be similar to 1Rj^ ph in the range under study, but higher at low 
redshifts, when dry mergers become more important. We can 
compare the gas-rich merger rate at z = 0.04 from the pre- 
vious fit, 2W0.04) = 0.52 x 10~ 4 Mpr 3 Gyr 1 , with the to- 
tal major merger rate inferred from the number of compani ons 
(N c = 0.011 ± 0.003) observed bv iDomingue et alJ (120091) at 
that redshift and stellar mass, % M = (0.8 + 0.4) x 10" 4 Mpr 3 
Gyr 1 (black square in Fig. [6j». The latter is higher than the for- 
mer and suggests that ~ 35% of the mergers at z = 0.04 are dry, 
in agreement with the expected tendency (see IWen et aP 2009 
for an estimation of the dry merger rate of luminous early-type 
galaxies at z ~ 0.1). 




Fig. 6. Major merger rate vs redshift for M* > 10 10 M remnant 
galaxies. White squares are from the present work by pair statis- 
tics, grey circles are f rom L09 by morpholog ical criteria, and the 
black square is from IDomingue et al.l (2009) by pair statistics. 
The black solid line is the best power-law fit to the z > 0.2 data 
and refers mainly to gas-rich major mergers (see text for details). 

Table 8. Major merger rate of M* > 10 10 M Q galaxies in 
GOODS-S 



Redshift 


<n mph 






(10~ 4 Mpcr 3 Gyr 1 ) 


(10- 4 Mpcr 3 Gyr 1 ) 


0.335 




1.1 ±0.7 


0.4 
0.725 


9 +2 6 


1.6 ±0.7 


0.725 


2 2 +L5 
-91 




0.975 


9 Q+1.6 
^-I 





6. The role of mergers in the evolution of 
intermediate-mass early types 

In a previous work, Lopez-Sani uan et al.l d2010l) have studied 
the number density evolution of early-type (ET, E/SO/Sa) galax- 
ies of log(M+/M ) > 10 at z < 1, finding that these galax- 
ies increase their comoving number density (pet) by a factor 
of 5 from z — 1 to the present. Comparing the evolution of 
the early-type population between z sup and zm < z sup , named 
P E e T w fonf'Zsu P ) = pET(zinf) - Pet(z sup ), against the gas-rich ma- 
jor merger rate in the same redshift range, we can define /et,m> 
the fraction of new early types that appear between z SU p and zm 
because of gas-rich major mergers: 



/ET.Mfonf , Zsup) - 



Piemfenf > Zsup) 
PETfonf ) - PET(Zsup) ' 



(29) 



where p^ m (zM,z S up) is the number density of gas-rich major 
merger remnants, 

Premfenf ' ^sup) — I 
Jin 



M - - 1 ' r m (o)(i+ zr 1 dz 



H () E{zY 



(30) 



where E(z) = -y^A + ^m(1 + z) 3 in a flat universe, and 51m(0) 
and n are the merger rate parameters in Eq. ( 1271 ). We as- 
sume that gas-rich merger remnants are early -type galaxies 
dNaab et al.ll2006al; iRothberg & Josephll2006allbt iHopkins et all 
120081 12009bh~ 
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Using the L09 gas-rich major merger rate, 
lLopez-Saniuan et al.1 (120101) infer that / et ,m(0, 1) = 17+'°% for 
early-type galaxies of log(M+/M ) > 10. As we have shown 
in Sect. [5] when we join the major merger rate by close pairs 
statistics from this work with the one in L09, derived from 
morphological criteria, the merger rate parameters change 
substantially from those derived by L09 using morphological 
information alone. However, the number density of major 
merger remnants remains similar in the two determinations, 
varying from p™ m (0, 1) = (8.5+f'jj) x 1(T 4 Mpc~ 3 with L09 
parameters to p™ m {0, 1) = (9.1+f*) x 1(T 4 Mpc~ 3 with the 
parameters from the present paper (Eq. 11281 ). which implies 
/et.m (0, 1) = 18*y 1 %. Both morphologies ( lLopez-Saniuan et al. 
2010) and close pairs (this work) therefore agree in yielding 
consistently low values for the fraction of early types that appear 
between redshifts 1 and due to major mergers. 



6.1. Estimating the minor merger rate 

Since gas-rich major mergers cannot expla in the rise in M+ > 
10 10 M o early-type galaxies since z ~ 1. lLo pez-Saniuan et alJ 
(2010) suggest minor mergers and secular processes (e.g., bars, 
disc instabilities, gas exhaustion, or morph ological quenching) 
as the main path i n that evolution (see also Bun dv et alj |2009b; 
lOesch et al. 2010). In Sect. l4.2l we find that the number of minor 
companions (N™) is roughly twice the number of major compan- 
ions (JV C M ). This implies that JV C (1/10 < yu < 1/3) = N™ - ~ 
N^ 1 , so we can estimate the minor merger rate (^R mm ), defined 
as the merger rate of galaxies with 1/10 < ft < 1/3. We denote 
= 2vm + Rimii = R x %m as the total (major + minor) merger 

rate. 

We apply Eq. (l26l i to obtain ^R mm . We assume that pa- 
rameters for minor companions are similar to those for major 
companions, except for the merg er time scale r pa j r . From N- 
body hydrodynamical simulations. lLotz et al.l d2009l) find T 



pair 



1.5 x T M for 5/T 

pair 



1 kpc < r p < 20h 1 kpc close pairs. With 
the previous assumptions we infer that % mm ~ 0.7 x %m> 
and ^R m ~ 1.7 x *&m- The latest cosmological models (e.g., 
Stewart et afl20 09a: Go nzalez-Garcfa et al.l2009l;lHopkins et al.1 
20T0F Dredict/? ~ 2, similar to our estimation R ~ 1.7. However, 



simulations refer to minor mergers selected by baryonic (gas + 
stellar) mass, while our estimation is from minor mergers se- 
lected by stellar mass. Although baryonic masses are more un- 
certain than stellar ones, it is interesting to compare the bary- 
onic and the stellar m erger rates, which could be quite different 
dStewart et al.ll2009bb . Moreover, the baryonic merger rate can 
be compared t o the fraction of ga laxies with distorted kinemat- 
ics at z ~ 0.6 dNeichel et al.ll2.Q08h . which are described well by 
merger simulations (//baryonic £1/5, Hammer et al. 20091). 



6.2. Mergers vs secular processes in the evolution of 
early-type galaxies 

Following the Sect. [6] steps with the combined major and mi- 
nor merger rates estimated in the previous section, we infer that 
~ 30%, and up to -50%, of the early-type (E/SO/Sa) galaxies of 
M+ > 10 10 M Q that appear since z ~ 1 may have undergone a 
major or minor merger event. This is an upper limit because a 
single minor m erger does not transfo rm a late-type galaxy into 
an early type (Hopkins et al. 2009b), but incr eases the Sersic 
index of the galaxy dEliche-Moral et aT] l2006t [Bournaud et al. 
12007b . This result suggests that the other ~ 50% of the new 



early t ypes appear due to secular processes (see also lBundv et al.l 
I2009bl) . 

That a large fraction of SO-Sb discs host pseudo-bulges 
dKormendv & Kennicutdl20"04T) also points to secular processes 
for the growth of the red sequence at z < 1. And the strong 
simil arity of disc and nuclear colours in di sc galaxies up to z ~ 
0.8 (iDommguez-Palmero & Balce lls 2008) likewise argues for a 
fading process that does not destroy the disc. The process leading 
to such an evolution may be suggested by the mentioned colour 
similarity of (the inner parts of the) disc and nucleus/bulge. Both 
star formation and the subsequent fading must be to some degree 
coordinated throughout the galaxy, so the qu enching of star for- 
mation may simply be du e to gas exhaustion (Zheng et al. 2007; 
iBauermeister et alj 2009). Interestingly, the star formation in the 
galaxy, leading to the growth of the central bulge component, 
may lead to this fading as the growth of the central potential sta- 
bilizes the gas compone nt preventing disc fragmentation (mor- 
phological quenching. iMartig et al.ll2009l) . Disc fragmentation 
contributes both to star formation and to the growth of the cen- 
tral b ulge, as shown by the evolution of chain and tadpole galax- 
ies (|Elmegreen et al.l2008l;ICeverin o et al. 2009; Bournau d et al.l 
2009). 

These results refer to intermediate-mass galaxies (M+ ~ 
4 x 10 10 M Q ) at z < 1, but the picture is different at higher stellar 
masses. On the one hand, M* > 10 11 M Q galaxies have higher 
pair fractions than less massive ones dde Ravel et al.l |2009; 
iBundv et al.l l2009a|). and re d pairs are more common at these 
masses dBundv et aTll2009al) . This suggests that dry mergers are 
an importan t process in the evolution of massive galaxies sinc e 
z ~ 1 (e.g.. iBell et ail 12001; iLin et al.ll2008l; lllbert et al.|[2010h 
On the other hand, the size ( Truiillo et al. 2007; Buitrago et al.l 
2008; Ivan der Wei et al.ll2008l: Ivan Dokkum et al.ll2010l) and ve- 
locity dispersion evolution (Cenarro & Truiilloll2009t) of M* > 
10" M early-type galaxies since z ~ 2 also supports the im- 
portance of mergers, especiall y the impact of minor mergers on 
the evolution of these system s dBezanson et al. 2009; Naab et al.l 
2009; Hopkin s et alJl2009ah . In addition, residual star formation 
in early-type (E/SO ) galaxies at z < 0.7 can also be explained 
by minor merging (Kav irai et al.ll2009 , 20101). This p roblem has 
also been analysed by Elic he-Moral et al.l (1201 Pali , who have 
modelled the evolution of luminosity function backwards in time 
for M+ > 10 11 M galaxies, selected according to their colours 
(red/blue/total) and their morphologies. They find that the ob- 
served luminosity function evolution can be explained naturally 
by the observed gas-rich and dry major merger rates and that 50- 
60% of today's E/S0 in this mass range were formed by major 
mergers at 0. 8 < z < 1, with a small number evolution since z = 
0.8 (see al so Crist obal-Hornillos et al.l 120091; Ivan der Wei et al.l 
2009; lllbert et al.l2010tlEliche-Moral et al.ll2010bl). The gas-rich 
major merger fractions assumed by Eliche-Moral et al. (2010a) 
are those from L09 for B-band selected galaxies (Mb < -20), 
which were obtained in a similar way than the morphological 
merger fractions used through pr esent paper. 

This makes M; ~ 10 11 M Q dPerez-Gonzalez et al. 1 120081) a 
transition mass at z < 1 : at higher masses mergers are an impor- 
tant process in the evolution of early-type galaxies, while other 
mechanisms dominate the observed evolution at lower masses. 



7. Conclusions 

In this paper we have developed a new method, which is based 
on the one used widely over spectroscopic surveys, to determine 
the mean number of companions per galaxy (A^ c ) over current 
spectro-photometric surveys. We tested our method in a local, 
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volume-limited sample from MGC spectroscopic catalogue. We 
find that the method provides reliable N c values when either 
photometric redshift errors are smaller than <x lnax or the spec- 

troscopic completeness of the sample is higher than /™", with 
these limits depending on the sample under study. For typical 
SDSS samples with tr Zptot = 0.02, we find f™£ ~ 0.2, while we 

find f™] ~ 0.4 for our GOODS-S catalogue. 

We studied the number of companions in Z?-band luminosity- 
selected samples, finding that 

- N c depends on search radius as expected from correlation 

function, N c oc ri, 3 , and their values are similar to those ex- 

p ' 

pected from a spectroscopic sample. 

- The values of 7V C for different luminosity selections are in ex- 
cellent agreement with those in the literature when the same 
selection criteria and pair definition are applied. We find that 
the number of companions decreases when luminosity in- 
crease and that the N c of Mb < -20 galaxies evolves faster 
(to = 2.5) than for Mb < -19 galaxies (to = 1.6). In addition, 
this evolution becomes slower, to = 0.7, when luminosity 
evolution is taken into account. 

We applied our new methodology to estimate the relation 
between the close major companions (N™, mass ratio ju > 1/3) 
and the close minor companions (/V™, mass ratio // > 1/10) of 
M+ > 10 10 M Q galaxies. We find that N™ is roughly two times 
N™. Studies in more extended fields are needed to better under- 
stand the minor-to-major ratio and their possible evolution with 
redshift. 

Finally, we compared the merger rates derived by close pairs 
with those by morphological criteria for M* > 10 10 M Q galaxies, 
finding that both are similar when the progenitor bias is taken 
into account. We estimate that the total (major + minor) merger 
rate is ~1.7 times the major merger rate. After comparing the 
merger rate with the observed structural evolution of GOODS- 
S galaxies, we infer that up to ~ 50% of the new early-type 
galaxies appeared since z ~ 1 could have undergone a merger 
event (major or minor), and we refer to star formation fading in 
the disc of spiral galaxies to explain the other ~ 50%. 

Our new methodology can be applied to extensive spectro- 
photometric surveys such as AEGIS, COSMOS, or VVDS, 
which have an order of magnitude more galaxies than the present 
catalogue. This will allow better determination of the number of 
major and minor companions and their dependence on redshift 
and other galactic properties, such as luminosity, mass, colour, 
or environment. 
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